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Abstract 

We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the 
intracluster medium in the merging galaxy cluster, 1E0657-56 (the Bullet cluster), using a series of N-body 
and hydrodynamic simulations. We find that the electron temperature at the shock layer associated with 
the X-ray sub peak (bullet) is quite different depending on the thermal relaxation model between electrons 
and ions; ~ 25 keV for the Coulomb thermal relaxation model and ^ 45 keV for the instantaneous thermal 
relaxation model in the simulations which reproduce the observed X-ray morphology. Furthermore, both 
of Fe XXV and Fe xxvi are overabundant compared with the ionization equilibrium state around the shock 
layer, and thus, the intensity ratio between Fexxv and Fexxvi Ka lines are significantly altered from 
that in the ionization equilibrium state. We also carry out the simulations with various sets of merger 
parameters, and discuss a possible range of the non-equilibrium effects in this system. Our results could 
be tested with future X-ray observations such as A.stro-H with better sensitivity in high energy band. 
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1. Introduction 

Cosmological shock heating induced by successive merg- 
ing of galaxies, galaxy groups, and clusters, are believed 
to be the main mechanism that have heated the X-ray 
emitting intracluster medium (ICM). Hydrodynamic sim- 
ulations of the structure formation have shown that shock 
waves in galaxy clusters has a Mach number of a few in the 
standard ACDM cosmology (e.g., Ryu et al. 2003), which 
heat the ICM at most several tens keV. Such very hot ICM 
has been found through X-ray emission and the Sunyaev- 
Zel'dovich (SZ) effect (e.g., RXC J1347.5-1144: Kitayama 
et al. 2004; Ota et al. 2008, Abell 3667: Nakazawa et 
al. 2009), suggesting that violent mergers are ongoing 
in the context of a standard scenario of the hierarchical 
structure formation in the universe. 

Cosmological shock waves are also important as labo- 
ratories of coUisionless plasma. For instance, turbulent- 
fiow motions induced via cascade of the vorticity gen- 
erated at the shocks amplify magnetic fields through 
the turbulent dynamo (e.g., Ryu et al. 2008). Related 
to the magnetic field, particle acceleration and nonther- 
mal emission are also interesting phenomena around the 
shocks (see Sarazin 2002 for a review). Distribution 
functions for shock heated electrons were recently stud- 
ied in the framework of the relativistic correction of 
the SZ effect (Prokhorov et al. 2011a; Prokhorov et 
al. 2011b; Prokhorov et al. 2011c). In last decade, 
a possible difference in temperature between electrons 



and ions, or the two-temperature structure, and the 
non-equilibrium ionization state of the ICM at the 
post-shock regions have been intensively investigated 
(Takizawa 1999; Takizawa 2000; Yoshida et al. 2005; 
Yoshikawa and Sasaki 2006; Cen, Fang 2006; Akahori 
and Yoshikawa 2008; Rudd, Nagai 2009; Akahori and 
Yoshikawa 2010; Wong et al. 2011). 

Akahori and Yoshikawa (2010) (hereafter AYIO) simu- 
lated merging galaxy clusters, relaxing the assumptions of 
both an ionization equilibrium and a thermal equiparti- 
tion between electrons and ions, and found that both as- 
sumptions are not justified around the shocks in merging 
galaxy clusters, because their relaxation timescales, order 
of 10^ yr for the electron density of 10~^ cm~'^, are not 
short enough compared with the merger timescale. They 
demonstrated that deviations from the ionization equi- 
librium and electron-ion equipartition states significantly 
affect an interpretation of the observational data such as 
the ICM temperature and metallicity. Therefore, studying 
the non-equilibrium effects is necessary not only to under- 
stand the relaxation mechanism in the plasma but also 
to correctly measure the properties of the shock-heated 
ICM with X-ray continuum emission of free electrons and 
X-ray line emissions of heavy elements. The effects of 
the non-equilibrium states in galaxy clusters have already 
been recognized with Suzaku X-ray observatory (Ota et 
al. 2008; Hoshino et al. 2010; Akamatsu et al. 2011). 

1E0657-56 or RXC J0658. 5-5556 {z = 0.296) is known 
as one of the hottest galaxy clusters. It hosts an X-ray 
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sub peak on ^ 500 kpc west side from a main peak in 
the ROSATBRl image (Tucker et al. 1998). Chandra ob- 
servations (Markevitch et al. 2002; Markcvitch 2006) have 
revealed that the X-ray sub peak has a "bullet" -like shape, 
that is why 1E0657-56 is called the Bullet cluster. Except 
for a region around the sub peak, average ICM tempera- 
ture was estimated to be kT = 17.4 ± 2.5 keV (Tucker et 
al. 1998) and 14.1 ±0.2 keV (Markevitch et al. 2002) from 
the fitting of the ASCA and Chandra data, respectively, 
which is relatively high compared with other observed 
clusters (e.g., Ota et al. 2006; Cavagnolo et al. 2009). 
Thermal SZ effect through such hot ICM was also ob- 
served (e.g., Zemcov et al. 2010). Optical-band obser- 
vations have shown that galaxy components are clearly 
offset from the associated X-ray peaks (e.g., Barrena et 
al. 2002), and the gravitational potential does not trace 
the distribution of the ICM but follows approximately 
the distribution of galaxies as expected for a collision- 
less dark matter component (Clowe et al. 2004; Clowe et 
al. 2006; Bradac et al. 2006). These features suggest that 
a violent merger of a sub cluster is ongoing in the 1E0657- 
56 system. A density jump at a bow-shock structure 
ahead of the bullet indicates a shock front with a Mach 
number of M = 3.0 ±0.4 (Markevitch 2006). Therefore, 
the 1E0657-56 system is an suitable target to investigate 
the non-equilibrium ionization state and the electron-ion 
two-temperature structure. Many theoretical works were 
devoted to numerical simulations of the 1E0657-56 sys- 
tem (Takizawa 2006; Milosavljevic et al. 2007; Springel & 
Farrar 2007; Nusser 2008; Mastropietro & Burkert 2008), 
and concluded that an ongoing merger of two galaxy clus- 
ters can reproduce the observed features such as very hot 
ICM, shocks, and displacement between X-ray and col- 
umn density peaks. These works, however, have assumed 
the ionization equilibrium of the ICM and the thermal 
equipartition between electrons and ions, and we are the 
first to address the non-equilibrium states of the ICM for 
further elucidation of the physical properties of the ICM 
in 1E0657-56 system. 

In this paper, we investigate an ionization state and 
temperature structure of the ICM in the 1E0657-56 sys- 
tem, and clarify to what extent deviations from ionization 
equilibrium and thermal equipartition between electrons 
and ions arise and affect an interpretation of the obser- 
vational data. We also discuss future prospects for the 
detectability of these non-equilibrium phenomena using 
the near future X-ray missions such as Astro- H. 

The rest of this paper is organized as follows. In sec- 
tion 2, we describe the model and setup of the simulations. 
The results are shown in section 3. We discuss parameter 
dependence of our results and effect of radiative cooling in 
section 4 and 5, respectively, and summarize our conclu- 
sions in section 6. In the case that the cosmological scaling 
is required, we assume the density parameter, flu = 0.24, 
the cosmological constant, Ha = 0.76, the baryon density 
parameter, f2b — 0.04, and the Hubble constant, Hq = 70 
km s~^ Mpc"-'^, unless otherwise specified. 



2. Numerical Simulation 

We carry out a set of N-body and SPH simulations 
of the 1E0657-56 system, relaxing the assumptions of 
the ionization equilibrium of the ICM and the thermal 
equipartition between electrons and ions. Detailed de- 
scriptions of our simulation code were already published 
in AYIO. Below, we briefly summarize our numerical 
schemes to treat the non-equilibrium ionization state and 
the thermal relaxation between electrons and ions as well 
as the setup of our simulations. 

2.1. Electron- Ion Two Temperature Structure 

Timescales on which each of electrons and ions achieves 
thermal relaxation through Coulomb scattering are much 
shorter than that between electrons and ions and a 
timescale of the merger in the cluster environment. 
Thus, we assume that electrons and ions always reach 
Maxwellian distributions with temperatures, Tg and Ti, 
respectively. Here, electrons and ions can have differ- 
ent temperatures just after experiencing shock heating, 
since most of the kinetic energy of the ICM is carried by 
ions at a shock front, and is preferentially converted to 
the thermal energy of ions in the post-shock regions (Fox, 
Loeb 1997). Afterward, thermal relaxation between elec- 
trons and ions proceeds toward the equipartition state. 

Physical processes responsible for the thermal relax- 
ation are still matter of debate. It is well known that 
the Coulomb scattering is not efficient enough to achieve 
electron-ion thermal equipartition in a timescale suffi- 
ciently shorter than typical merger timescale, ~ 10^ yr, 
and it is argued that plasma waves are able to attain 
the thermal relaxation between electrons and ions more 
quickly than the Coulomb scattering. In this work, in 
order to bracket the plausible range of theoretical un- 
certainties regarding the thermal relaxation processes be- 
tween electrons and ions, we present the results with two 
different thermal relaxation models: "single-temperature 
runs" in which the equipartition of thermal energy be- 
tween electrons and ions is achieve instantaneously, and 
"two-temperature runs" which adopt only the Coulomb 
scattering as a physical process for thermal relaxation be- 
tween electrons and ions. 

In the two-temperature runs, we solve an ordinary en- 
ergy equation for mixed fluid with a mean temperature, 
T = (ucTc + n-iTi)/ [ric + rij), where and rij the number 
density of electrons and ions, respectively, and another 
equation for the time evolution of electrons which can be 
reduced to 

dt tci u 

where Te = Te/T and T{=T{/T are the dimensionless tem- 
peratures of electrons and ions, respectively, normalized 
by the mean temperature, t^i the Coulomb thermal relax- 
ation timescale between electrons and ions (equation (3) 
of AYIO), and u and Qsh are the specific thermal energy 
and the shock heating rate per unit mass, respectively. 
We solve equation (1) in the same manner described in 
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Takizawa (1999). 

The effect of radiative cooling of the ICM is con- 
sidered in our simulations, in which only the thermal 
bremsstrahlung emission is considered since the ICM tem- 
perature is always high enough a few keV). We com- 
pute the cooling rate by using T in the single-temperature 
runs and Tg in the two-temperature runs to incorporate 
the effect of two-temperature structure of the ICM con- 
sistently. 

2.2. Non-Equilibrium Ionization State 

The Lagrangian time evolution of ionization fractions of 
heavy elements, C, N, O, Ne, Mg, Si, S, and Fe, as well as 
hydrogen and helium is computed for each SPH particle 
by solving rate equations: 

fc=l i=i+l 

where j is the index of a particular ionization stage con- 
sidered, Z the atomic number, fj the ionization fraction 
of an ion j, Si^j the ionization rate of an ion j with the 
ejection of i electrons, and aj is the recombination rate 
of an ion j. The ionization processes include collisional. 
Auger, charge-transfer, and photo-ionizations, and recom- 
bination processes are composed of radiative and dielec- 
tronic recombinations. We solve equation (2) in the same 
manner described in Yoshikawa and Sasaki (2006). 

Ionization and recombination rates are calculated by 
utilizing the SPEX ver 1.10 software package^. In the two- 
temperature runs, the reaction rates are computed using 
the electron temperature, Te, calculated from equation (1) 
in order to incorporate the effect of the two-temperature 
structure. As for the single-temperature runs, the reac- 
tion rates are computed using the mean temperature, T . 
All ionization states of the atoms denoted above are con- 
sidered in calculations of X-ray surface brightness. But in 
the following we only focus on the ionization state of iron, 
because most of ions except iron are fully-ionized for the 
ICM with the temperature higher than several keV in the 
1E0657-56 system. 

2.3. Simulation Setup 

As an initial condition of the 1E0657-56 system, two 
spherically symmetric galaxy clusters are set up as fol- 
lows. The virial radius, r200j within which the mean 
mass density is 200 times the present cosmic critical den- 
sity, and the virial mass, M200, enclosed within r2oo of 
a main (massive) cluster are set to r2oo = 2.36 Mpc and 
-^200 = 1.5 X 10^^ A/q, respectively, and those of a sub 
(less massive) cluster are 1.21 Mpc and 2.5 x 10^^ Mq, re- 
spectively. We adopt the hydrostatic model (AYIO) with 
the NFW density profile (Navarro et al. 1997) for a dark 
matter halo, and the /3-model density profile (Cavaliere 
& Fusco-Femiano 1976) with ;9 = 2/3 for an ICM com- 
ponent, and both of dark matter and SPH particles are 
initially distributed out to r2oo- The scale radii, rg, of the 

^ ( http://www.sron.nl/divisions/hea/spcx/ ). 



NFW profile arc set so that the concentration parameters 
c = rs/r2oo 5 and 8 for the main and sub clusters, re- 
spectively, and the core radius of the /3-modcl is set to 
fc —T&l'2. for each cluster. 

The main cluster is composed of lO'' dark matter parti- 
cles and the same number of SPH particles, while those of 
the sub clusters are proportional to their mass ratio. The 
corresponding spatial resolution of SPH is 12 kpc at the 
ICM density of 10"^ cm"'"' (at the bullet) and 26 kpc at 
10~^ cm~^ (at the shock front). 

Our nominal run that reproduces the gross observed 
features of the 1E0657-56 system is organized as fol- 
lows. Initially, two galaxy clusters contact each other 
at their outer edge, ^200- Their initial relative veloc- 
ity is set to 3000 km/s, and the impact parameter is 
h = 0.236 Mpc. Such an initial relative velocity and a non- 
zero impact parameter are required to reproduce a shock 
front with a Mach number of M ~ 3 and asymmetric X- 
ray surface brightness (section 4; see also Mastropietro & 
Burkert 2008). As for the metallicity, a spatially uniform 
metallicity of 0.2 times solar abundance is assumed. In 
addition to that, we also assume Tc = Tj = 1 and an ion- 
ization equilibrium state at the start of the simulations. 

Compared with the previous hydrodynamic simulations 
(Springel & Farrar 2007; Mastropietro & Burkert 2008), 
the numbers of SPH and dark matter particles are ■^S-IO 
times larger, so that the spatial resolution is roughly twice 
better. While Mastropietro & Burkert (2008) adopted the 
larger concentration parameters and smaller main cluster 
mass, we adopt the larger concentration parameters but 
the same main cluster mass as those adopted in Springel 
& Farrar (2007). This is because our parameters better 
reproduce the observed temperature profile at the pre- 
shock region in front of the bullet. While Springel & 
Farrar (2007) neglected the effect of radiative cooling, we 
demonstrate that considering radiative cooling is impor- 
tant since the radiative cooling significantly affects the 
width between the shock front and the contact disconti- 
nuity associated with "the bullet" (section 5). Although 
our nominal run reproduces an overall structure of the 
1E0657-56 system well, as a practical interest, we also 
carried out the simulations with other sets of parameters, 
and discuss a possible range of the non-equilibrium effects 
(section 4). 

3. Result 

In the rest of this paper, we present the result after 
1.12 Gyr has elapsed from the initial condition. At this 
epoch, the sub cluster already penetrates the center of 
the main cluster, and the separation between the centers 
of the main and sub clusters reaches ^ 0.72 Mpc, which is 
consistent with that of the 1E0657-56 system (e.g., Clowe 
et al. 2004), assuming that the collision plane is perpen- 
dicular to the line-of-sight (the linc-of-sight velocity differ- 
ence of galaxies is estimated to be only 600 km/s; Barrena 
et al. 2002). 
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Fig. 1. X-ray surface brightness in the 0.8—4.0 keV band. 
The projected mass density (dark matter -|- ICM) is also over- 
laid (white contours). The core of the sub cluster (the right 
cluster) is moving from left to right hand side and has already 
penetrated the main cluster (the loft cluster). 

3.1. Overall Structure of Simulated 1E0657-56 

We first present the overall structure of the simulated 
clusters in order to make sure that our simulation repro- 
duces the gross observed features of the 1E0657-56 sys- 
tem, and to locate the regions where the effects of the 
non-equilibrium ionization state and the two-temperature 
structure remarkably take place. 

Figure 1 depicts the 0.8-4.0 keV band X-ray surface 
brightness map and the contours of the projected mass 
density (dark matter -|- ICM) integrated along the line- 
of-sight. Our nominal run nicely reproduces the over- 
all structure of the 1E0657-56; for instance, distinct 
two X-ray peaks, positions and shapes of shock front 
and contact discontinuity, and displacement between X- 
ray and projected mass density peaks associated with 
the sub cluster, are consistent with the observations 
(Markevitch 2006; Clowe et al. 2004). 

Figure 2 shows the temperature maps on the collision 
plane of the two galaxy clusters. There exists very hot 
ICM heated from ~ 8 keV to more than ^ 30 keV by a 
shock wave with a Mach number of M 2-3 in front (the 
right hand side) of the cold core (the bullet) of the sub 
cluster. The shock layer extends toward the outskirts of 
the clusters, which has a Mach number oi M ^ 1.5-2.5 and 
heats the ICM to ^ 15-20 keV. In the two-temperature 
run, we can see that the electron temperature is typi- 
cally ~ 30-50 % lower than the mean temperature be- 
hind the shock (the bottom panel of figure 2). The two- 
temperature structure can not be seen inside the bullet, 
because the timescale of thermal relaxation between elec- 
trons and ions are much shorter than the merger timescale 
due to the high density and low temperature. We have an- 
other high temperature region behind the bullet heated by 
the compressive flows. 

Figures 3 and 4 show the ratios of Fe xxv and Fe xxvi 
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Fig. 2. The mean temperature of ICM, T, (top) and the ratio 
of the electron temperature relative to the mean temperature, 
Tc/T, (bottom), on the collision plane of the two galaxy clus- 
ters. Contours show the Mach number on the collision plane 
from 1.5 by 0.5, and the innermost contour represents the 
Mach number of 2.5, estimated by the shock-capture method 
described in Pfrommer et al. (2006). A white horizontal line 
indicates the location of the one-dimensional profile discussed 
in figure 7. 

fractions relative to those in the ionization equilibrium 
state, respectively. Here, in the two-temperature run, 
ionization fractions in the ionization equilibrium state are 
computed using the mean temperature. We find that both 
Fexxv and Fexxvi are overabundant by 1.2-3.0 behind 
the shock layer, and the deviations from the equilibrium 
value are larger in the two-temperature run. The over- 
population of Fe XXV and Fe xxvi fractions can be under- 
stood as follows. The shock heats the ICM in the post- 
shock region from several keV to a few tens keV, and both 
of Fe xxv and Fe xxvi fractions in the post-shock region 
decrease with time toward the equilibrium value because 
their fractions are maximized at ^ 3 keV and ~ 10 keV, re- 
spectively, in the ionization equilibrium state (see figure 4 
of AYIO). However, the ionization of Fexxv and Fexxvi 
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Fig. 3. The ratio of the ionization fraction of Fo XXV relative 
to that in the ionization equilibrium state on the collision 
plane of the two galaxy clusters. Top and bottom panels 
show the results of the single- and two-temperature runs, re- 
spectively. Black contours are the same as figure 2. 

toward higher-ionized levels is not quick enough to catch 
up with the ionization equilibrium state. Therefore, these 
fractions are higher than that in the ionization equilib- 
rium state. Since the ionization rates of ions to the higher 
ionization levels depend on the electron temperature, and 
the rates are significantly smaller in the two-temperature 
run, a delay of the heating of electrons at the shock leaves 
these fractions higher than that in the single-temperature 
run. 

To see an observational appearance of the result of the 
nominal runs, we calculate the spectroscopic-like temper- 
ature (Tsi, Mazzotta et al. 2004) and the ratio, R/Rcq, 
as an observational tracer for the non-equilibrium ioniza- 
tion state (AYIO), where R is the X-ray intensity ratio 
between rest-frame 6.6-6.7 keV and 6.9-7.0 keV energy 
bands in which the Fe xxv Ka and Fe xxvi Ka emissions 
are predominant, respectively; 




-1 

-0.5 0.5 1 

X [Mpc] 



Fig. 4. Same as figure 3 but for the ionization fraction of 
Fexxvi. 

7(6.6- 6.7 keV) 

7(6.9- 7.0 keV)' ^ ' 

and 7?oq is the intensity ratio defined above but in the 
ionization equilibrium state. It should be noted that the 
quantity, 7?/7?oq, is independent on the ICM metallicity. 
The results are shown in figures 5 and 6. In terms of 
the spectroscopic-like temperature, the temperature of the 
ICM at the shock in front of the bullet is roughly ^15 keV 
in the two-temperature run, and lower than in the single- 
temperature run. R/Rcq departs from unity (^ 1.1-1.2) 
both in the single- and two temperature runs, though its 
difference between the two runs is not significant. In both 
single- and two-temperature runs, the deviation of the in- 
tensity ratio is more significant in the outskirts than in the 
central part of the clusters, since the timescale to reach the 
ionization equilibrium state is longer in the lower-density 
regions. 
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Fig. 5. The spectroscopic-like temperature (Mazzotta et al. 
2004). Top and bottom panels show the results of the single- 
and two-temperature runs, respectively. 



Fig. 6. R/Rcq (see equation (3)). Top and bottom panels 
show the results of the single- and two-temperature runs, re- 
spectively. 



3.2. Detailed Structures of Bullet and Shock 

The shock front ahead of the bullet is suited to the 
study of the non-equilibrium ionization state and the two- 
tcmpcrature structure of the ICM, since the strong shock 
with a Mach number of 3.0 ± 0.4 is reported (Markevitch 
2006) and the observational signatures of the departure 
from the ionization equilibrium and the two-temperature 
structure are expected to be relatively strong. In this 
subsection, we present the detailed analyses of the shock 
ahead of the bullet, and how such observational signatures 
can be detected in X-ray observations. 

Figure 7 shows the one dimensional profiles of the 
ICM across the bullet and the shock front on a line 
on the collision plane, the location of which is indi- 
cated by a white line in figure 2. We can see that the 
shock velocity, the relative velocity of the ICM between 
pre- and post-shock regions, is about 5000 km/s (fig- 
ure 7(b)), which is consistent with the previous simula- 
tions (Takizawa 2006; Milosavljevic et al. 2007; Springel 



& Farrar 2007; Mastropietro & Burkert 2008), The 
Mach number of the shock is '^2.8, and consistent with 
the estimate based on X-ray observations by Chandra 
(Markevitch 2006). 

In the two-temperature run, we find that the electron 
temperature between the shock front and the contact 
discontinuity is Tq ~ 25 keV and much lower than the 
mean temperature, f ^ 45 keV. In the single-temperature 
run, the temperature of the ICM is almost the same as 
the mean temperature in the two-temperature run, and 
the profiles of the ICM density and velocity are also 
unchanged regardless of the thermal relaxation model 
adopted. This means that the difference in the cooling 
rate due to the difference in the electron temperature does 
not affect the overall hydrodynamic properties of the ICM. 

We can see that the ionization state of iron strongly de- 
viates from the ionization equilibrium state between the 
shock front and the contact discontinuity, and the ratios 
of Fexxv and Fexxvi fractions relative to those in the 
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Fig. 7. One dimensional profiles across the shock front on 
the collision plane. Thin and thick lines show the results of 
the single- and two- temperature runs, respectively. Panels 
from the top to the bottom show (a) the ICM density, (b) 
the ICM velocity along the axis, (c) the average temperature 
(solid line) and the electron temperatures in the two-temper- 
ature run (dashed line), (d) the ratio of ionization fractions 
of Fe XXV relative to that in the ionization equilibrium state, 
and (e) the ratio for Fexxvi. Two vertical lines indicate the 
location of the shock front (right) and the contact disconti- 
nuity (left). The coordinate of the horizonal axis shows the 
distance measured from the shock front. 

ionization equilibrium state reach 5 and 2.5, respec- 
tively (figures 7(d) and 7(e)). The amount of the depar- 
ture from the ionization equilibrium state depends on the 
adopted thermal relaxation model and larger in the two- 
temperature run. 

Figure 8 shows the profiles of X-ray surface bright- 
ness in the 0.8-4.0 keV band (panel (a)), spcctroscopic- 
like temperature (panel (b)), and R/Rcq (panel (c)) pro- 
jected along the line-of-sight at the same region as figure 7. 
In these projected quantities, the signatures of the two- 
temperature structure and the non-equilibrium ionization 
state are contaminated by the contributions of the fore- 
ground and background ICM. Nevertheless, we can see 
such signatures in figure 8 to some extent. In the two- 
temperature run, the maximum spectroscopic-like tem- 
perature between the shock front and the contact discon- 
tinuity is ^ 20 keV and significantly lower than that in 
the single-temperature run (^ 25 keV). R/Rcq also devi- 
ates from unity at the same region, indicating the non- 
equilibrium ionization state in this region, though the dif- 
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Fig. 8. Same as figure 7 except that panels from the top to 
the bottom show (a) the surface brightness in the 0.8-4.0 keV 
energy band, (b) the spectroscopic-likc temperature, and (c) 
R/Rcq projected along the line-of-sight, respectively. 

ference between the single- and two-temperature runs is 
not significant. 

It should be noted that, due to the poor capability of 
the SPH scheme to resolve the shock front, the discon- 
tinuities of physical quantities such as density and tem- 
perature across the shock front arc numerically smeared 
in our simulations, and that the actual jumps of density, 
velocity and temperature of the ICM should be sharper in 
"reality". The shaded regions in figures 7 and 8 indicate 
the locations where such numerical smearing is effective. 
If the jumps of physical quantities at the shock are ideally 
resolved, it is expected that the difference in tempera- 
ture between electrons and ions at the post-shock regions 
and the departure from the ionization equilibrium state 
are larger at the post-shock region. Therefore, the differ- 
ence of projected profiles of the spectroscopic-like temper- 
ature between the single- and two-temperature runs would 
be larger, and thus the difference of R/Rcq between the 
single- and two-temperature runs would be also more sig- 
nificant. 
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4. Dependence on the Initial Condition 

In this section, in order to figure out parameter depen- 
dence of our results, we compare various merger simula- 
tions with different sets of parameters. 

First, we briefly summarize morphological differences of 
the simulated clusters. The foUowings are notable depen- 
dences on the initial parameters that we confirmed: 

• An overall structure of the observed features such as 
the bullet, the shock layer, and the peaks of X-ray 
surface brightness and projected mass density are 
best reproduced with the concentration parameter, 
c = 5, of the main cluster. The runs with c = 3 (e.g., 
Springel & Farrar 2007) give similar results, but pro- 
duce X-ray and projected mass density maps of the 
main cluster slightly more elongated along north- 
south direction. Difference in the simulations with 
c = 5 and those with c = 6 adopted by Mastropietro 
& Burkert (2008) is almost negligible. 

• The displacement between the peaks of X-ray sur- 
face brightness and projected mass density in the 
sub cluster (the bullet) are best reproduced with 
the concentration parameter, c = 8, of the sub clus- 
ter (Mastropietro & Burkert 2008). For c> 8, the 
bullet gets close to the projected density peak, and 
the displacement becomes smaller. 

• As pointed out by Mastropietro & Burkert (2008), a 
non-zero impact parameter is necessary to reproduce 
the asymmetric X-ray surface brightness. However, 
the runs in which the sub-cluster gets through the 
south side of the main cluster can not simultane- 
ously reproduce the observed X-ray surface bright- 
ness and projected mass density. Therefore, we 
adopt the sub-cluster's orbit which passes through 
the north side of the main cluster. 

• The initial relative velocity of 3000 km/s is required 
to reproduce the shock layer with a Mach number 
of M 3 in front of the bullet. Similar requirement 
has been also indicated in other numerical studies 
on the bullet cluster (Takizawa 2006; Milosavljevic 
et al. 2007; Springel & Farrar 2007; Mastropietro & 
Burkert 2008). According to Hayashi, White (2006), 
an encounter with the virial velocity of the main 
cluster, V200 ^ 1700 km/s, and the relative veloc- 
ity, Vsuh ~ 3000 km/s, is probable with a rela- 
tive likelihood of f ^ 0.03, based on the cosmolog- 
ical structure formation simulations (see equation 
(1) of Hayashi & White 2006, where we adopted 
'I'lOpercent = 1.55, and a — 3.3). However, it should 
be noted that the probability for such a high ini- 
tial relative velocity is still controversial: Lee, 
Komatsu (2010) recently claimed that the initial 
relative velocity of 3000 km/s at i?200 is very un- 
usual with the prediction of the concordance ACDM 
model. Even if wc admit a M ^ 2.5 shock as the 
lower limit of the observed Mach number, the re- 
quired initial relative velocity is ~2500 km/s, and is 
still peculiarly high. 
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• The contrast of the peak brightness relative to that 
of the ambient gas is best reproduced with the ini- 
tial mass ratio of 6 : 1 between the main and the sub 
clusters. With decreasing the mass ratio, the X-ray 
peak of the main cluster becomes dimmer (see e.g., 
3 : 1 runs of Mastropietro & Burkert 2008). Note 
that the mass ratio also affects the Mach number of 
the shock; the run with the mass ratio of 4 : 1 pro- 
duces the shock with the Mach number larger than 
3 for the initial relative velocity of 3000 km/s, and 
allows us to adopt smaller initial relative velocity to 
reproduce the M 3 shock. 

Using the various merger simulations with different sets 
of parameters, let us evaluate a possible range of the non- 
equilibrium effects in the 1E0657-56 system. Based on 
the simulations in this and previous works (Springel & 
Farrar 2007, Mastropietro & Burkert 2008), we consider 
the results with the initial mass ratio, 4:1, 6:1, or 10 : 1, 
and the concentration parameter c = 3 or 5 for the main 
cluster and 8 or 10 for the sub cluster, which are con- 
sistent with the observational features of the 1E0657-56 
system to some extent. To identify the snapshot that cor- 
responds to the "current" epoch in the time sequences of 
the simulations, we adopt the separation of two peaks of 
the projected mass density. The other parameters such as 
the mass of the main cluster and the impact parameter 
are the same as our nominal run except that we use one 
eighth of the number of particles (the smoothing length 
is thus twice as large). We adopt the same initial relative 
velocity of 3000 km/s regardless of the mass ratio, since, 
in our prescription to determine the initial relative veloc- 
ity described in Sarazin (2002), it is proportional to the 
square root of the mass ratio and has weak dependence 
on the mass ratio. 

In table 1, we present the physical properties of the 
shock front associated with the bullet for the runs with 
the simulation parameters different from our nominal run. 
We find that R/R^q = 1.1-1.4 in the two-temperature 
runs, with the Mach number ranging from M = 2.5 to 
3.4. We also find that T = 33-50 kcV and = 9-24 keV, 
or Tc/T = 0.2-0.7. In other words, the observed electron 
temperature at the shock lower than ^ 25 keV suggests 
that the Coulomb collision is a predominant process of 
thermal relaxation between electrons and ions and that 
there exists the two-temperature structure at the shock 
front. On the other hand, the observed electron tem- 
perature higher than ~ 30 keV prefer the instantaneous 
thermal relaxation at the shock front. 

The Chandra observation suggests that the deprojected 
electron temperature behind the shock is over ^ 30 keV 
(Markevitch 2006), which is higher than the constraint 
for the existence of the two-temperature structure we ob- 
tained in the two-temperature runs. This suggests that, 
although there are still uncertainties due to the spatial res- 
olution in both observations and simulations, the instan- 
taneous thermal relaxation model is preferable for repro- 
ducing the observed temperature structure at the shock. 
Markevitch (2006) discussed that the simple model of 
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Table 1. Physical properties of the shock in front of the bul- 
let for the runs ^^cXcYmZ" where X and Y denotes the con- 
centration parameters of the main and sub clusters, respec- 
tively, and the mass ratio of the two clusters is Z : 1. f and 
To are in units of keV. 



run ID 


M 


T 


Tc 


Tc/T 


R/ Rcq 


c03c08m04 


3 





44 


21 





47 


1 


25 


c03c08m06 


3 





43 


21 





50 


-1 
1 


ZD 


c03c08ml0 


2 


6 


33 


20 





62 


1 


24 


c03cl0m04 


3 


1 


49 


20 





40 


1 


29 


c03cl0m06 


3 


4 


41 


13 





31 


1 


36 


c03cl0ml0 


3 


2 


34 


8.3 





26 


1 


19 


c05c08m04 


2 


9 


47 


24 





51 


1 


27 


c05c08m06 


3 





43 


23 





54 


1 


19 


c05c08ml0 


2 


5 


34 


22 





65 


1 


18 


c05cl0m04 


3 


1 


50 


22 





43 


1 


23 


c05cl0m06 


3 





49 


16 





34 


1 


14 


cOSclOmlO 


3 


2 


42 


9.9 





23 


1 


16 



electron-ion thermal relaxation only through the Coulomb 
scattering is inconsistent with the observed data and is ex- 
cluded at a 95% confidence level, implying that the ther- 
mal equilibration between electrons and ions should be 
much faster than that by the Coulomb scattering. We, 
however, should notice that the accuracy of the observed 
electron temperature by Chandra is somewhat limited in 
the high energy band more than ^ 10 keV. Thus, this is- 
sue should be addressed by the observations with future 
X-ray satellites such as Astro-H with better sensitivity in 
high energy band. Another plausible way to measure the 
temperature of very hot ICM should be a joint analysis 
with the X-ray and SZ data (Kitayama et al. 2004). The 
analysis of the Bullet cluster that is based on the high- 
resolution X-ray and SZ observations such as Chandra 
and Herschel-SPIRE will provide us new constraints on 
very hot thermal and/or non-thermal particle populations 
(Prokhorov, Akahori, Yoshikawa, et al., in preparation). 

5. EfTect of Radiative Cooling 

Finally, we discuss the effect of radiative cooling on the 
properties of the simulated bullet clusters. 

As partly discussed in section 3.2, the thermal relax- 
ation model docs not significantly change the overall hy- 
drodynamic properties of the ICM. This is because the 
radiative cooling is ineflScient in the pre-shock region and 
the shock layer, where radiative cooling timescale is longer 
than the age of the universe, regardless of the difference in 
the cooling rate due to the difference in the electron tem- 
perature. As for the bullet, radiative cooling is efficient, 
but there arc no two-temperature structure due to a short 
relaxation timescale between electrons and ions, so that 
the effect of radiative cooling is identical there. Therefore, 
morphology of the X-ray surface brightness does not de- 
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Fig. 9. One dimensional profiles across the bullet and the 
shock on the collision plane for the two-temperature runs. 
Thin and thick lines show the results of the non-cooling and 
cooling runs, respectively. Panels from the top to the bottom 
shows (a) the ICM velocity along the axis, (b) the ICM den- 
sity, and (c) the average temperature (solid) and the electron 
temperature (dashed). 

pend significantly on the thermal relaxation model, and 
thus we basically obtain the same best fit parameters for 
different thermal relaxation model. 

Considering the effect of radiative cooling is, however, 
quite important to determine the location of the bullet in 
the 1E0657-56 system. In order to demonstrate the effect, 
we carried out the control run in which we neglect the ef- 
fect of radiative cooling, and compared it with the cooling 
run that we have shown. The results after 1.12 Gyr has 
elapsed for the two-temperature runs are shown in figure 
9, where the axes of the profile for two runs are the same 
as that shown in figure 2. We can see a delay of the prop- 
agation of the bullet and the shock in the non-cooling run. 
Moreover, the width between the shock front and the con- 
tact discontinuity becomes wider in the non-cooling run 
than that in the cooling run. This may be ascribed to 
the fact that the gravitational binding force is weaker, i.e. 
ram pressure becomes more significant, in the non-cooling 
run, since the bullet less concentrates compared with the 
cooling run (by a factor of ^ 1.6 at the present). Our re- 
sult suggests that considering the effect of radiative cool- 
ing is important to correctly constrain the initial relative 
velocity from the measurement of the width between the 
shock front and the contact discontinuity (Mastropietro 
& Burkcrt 2008). 

6. Summary and Conclusion 

We conduct a set of A^-body/SPH simulations of merg- 
ing galaxy clusters which reproduce observed X-ray sur- 



10 



T. Akahori and K. Yoshikawa 



face brightness and projected mass density distribution of 
the 1E0657-56 system to investigate the non-equilibrium 
ionization state and clcctron-ion two-tcmpcrature struc- 
ture in the 1E0657-56 system, and to assess their de- 
tectability in future X-ray missions. Different from previ- 
ous numerical works of the 1E0657-56 system, we, for the 
first time, relax the both assumptions of ionization equi- 
librium and thermal equipartition between electrons and 
ions. 

Our nominal run, which reproduces the various ob- 
served features of the 1E0657-56 system fairly well, is an 
merger of two galaxy clusters with a mass ratio of 6:1 and 
an initial relative velocity of 3000 km/s. It is found that 
departure from the ionization equilibrium is significant at 
the shock front with a Mach number of M ^ 3 associated 
with the penetrating core of the sub cluster. When con- 
sidering only the Coulomb scattering as a physical process 
for the thermal relaxation between electrons and ions, we 
find that the two-temperature structure, or the difference 
in temperature between electrons and ions, are also sig- 
nificant at the shock front. 

Comparisons of the X-ray and column density profiles 
with those in the observations gave possible range of sim- 
ulation parameters. Within the range, at the shock asso- 
ciated with the bullet we found that R/Rcq (equation (3)) 
clearly deviates from unity and the electron temperature 
is much lower than the mean temperature of the ICM in 
the two-temperature run. Effect of radiative cooling is 
important to determine the location of the bullet, where 
the overall hydrodynamic properties of the ICM is not 
significantly altered by the different efhciency of radiative 
cooling due to the different temperature structure. 

It is found that, in the two-temperature run, the 
spectroscopic-like temperature projected along the line- 
of-sight between the shock front and the contact dis- 
continuity is significantly lower than that in the single- 
temperature run, which suggests that precise measure- 
ments of the maximum electron temperature in this re- 
gion could provide strong constraints on thermal relax- 
ation processes between electrons and ions. Such observa- 
tions can be available with a high sensitivity X-ray detec- 
tor in high energy band and also with a good spatial reso- 
lution to resolve the relatively tiny area between the shock 
front and the contact discontinuity. Observations of ther- 
mal Sunyaev-Zel'dovich effect with relatively good spa- 
tial resolution can provide an independent measurement 
of electron temperature in this region. A good spectro- 
scopic resolution in X-ray observations is also important 
in studying the two-temperature structure of the ICM, be- 
cause the ion temperature is only estimated from the de- 
tailed line profiles of emission lines, and would be achieved 
by X-ray calorimeters on board the forthcoming satellites. 
Our results also suggest that the current X-ray observa- 
tions can lead to a biased result on the estimation of the 
metallicity if one assumes that the ICM is in the ioniza- 
tion equilibrium state, because the separate detection of 
emission lines of Fe xxv and Fe xxvi is still difficult with 
the current X-ray observational facilities due to a lack of 
the spectroscopic resolution of the CCDs. Therefore, a 



good spectroscopic resolution is also needed to correctly 
measure the metallicity in the shock heated ICM. 
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